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ABSTRACT 



High resolution far ultraviolet spectra of the early-type stars 7 2 Vel and £ Pup were 
obtained to measure the interstellar deuterium abundances in these directions. The 
observations were made with the Interstellar Medium Absorption Profile Spectrograph 
(IMAPS) during the ORFEUS-SPAS II mission in 1996. IMAPS spectra cover the 
wavelength range 930-1150A with A/AA ~ 80,000. The interstellar D I features are 



resolved and cleanly separated from interstellar H I in the Ly<5 and Lye profiles of 



both sight lines, and also in the Ly7 profile of £Pup. The D I profiles were modeled 



using a velocity template derived from several N I lines in the IMAPS spectra recorded 



at higher signal-to-noise. To find the best D I column density, we minimized \ f° r 



model D I profiles that included not only the N(D I) as a free parameter, but also 
Q H i the effects of several potential sources of systematic error which were allowed to vary 

as free parameters. H I column densities were measured by analyzing Lya absorption 
+3 . profiles in a large number of IUE high dispersion spectra for each of these stars 

and applying this same x 2 -minimization technique. Ultimately we found that D/H 
= 2.18±gjf x 10~ 5 for 7 2 Vel and 1.42±°;|| x 10~ 5 for CPup, values that contrast 
markedly with D /H derived in Paper I for 5 Ori A (the stated errors are 90% confidence 



limits). Evidently, the atomic D/H ratio in the ISM, averaged over path lengths 
of 250 to 500 pc, exhibits significant spatial variability. Furthermore, the observed 
spatial variations in D/H do not appear to be anticorrelated with N/H, one measure 
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of heavy element abundances. We briefly discuss some hypotheses to explain the D/H 
spatial variability. Within the framework of standard Big Bang Nucleosynthesis, the 
large value of D /H found toward 7 2 Vel is equivalent to a cosmic baryon density of 
Q B h 2 = 0.023 ± 0.002, which we regard as an upper limit since there is no correction 
for the destruction of deuterium in stars. 

Subject headings: ISM: Abundances — Cosmology: Observations — ISM: Evolution 
- Ultraviolet: ISM — Stars: Individual (7 2 Vel, ( Pup) 

1. INTRODUCTION 

The abundance ratio of atomic deuterium to hydrogen (D/H) in interstellar gas is widely 
regarded as an important tracer of galactic chemical evolution (Audouze & Tinsley 1974; Boesgaard 
& Steigman 1985; Tosi et al. 1998) and a key discriminant of the cosmic baryon-to-photon 
ratio (if) in Big Bang Nucleosynthesis (BBN, Walker et al. 1991). A standard interpretation 
is that D should not be produced in significant quantities in astrophysical sites other than the 
Big Bang (Epstein, Latimer, & Schramm 1976). The generally accepted viewpoint is that while 
D is produced primordially, its destruction takes place when some of the gas is cycled through 
stars. The uncertainties surrounding this process represent a stumbling block in arriving at the 
primordial value. For this reason, a measurement of D/H is often regarded as a lower limit to the 
primordial ratio, and this may be translated into an upper limit to rj (for BBN, larger D/H implies 
lower rj). 

Rogerson & York (1973) made the first measurement of the atomic D/H abundance ratio 
in the interstellar medium (ISM) on the line of sight toward (3 Cen. Measurements with the 
Copernicus satellite toward an additional 14 stars (100 < d < 500 pc) found that D/H values 
cluster around 1.5 x 10~ 5 , but with a dispersion that, in some cases, seemed to exceed the stated 
uncertainties (see review by Vidal-Madjar & Gry 1984). While a simple interpretation of the 
Copernicus data suggested that the D/H measurements revealed spatial variations, the reality of 
these differences has been difficult to substantiate due to concerns stemming from the somewhat 
inadequate resolution of Copernicus (15 km s~ x FWHM) for this purpose. 

Using the Goddard High Resolution Spectrograph (GHRS) and Space Telescope Imaging 
Spectrograph (STIS) on the Hubble Space Telescope (HST), several measurements of D/H in the 
local interstellar medium (LISM, d < 100 pc) have been made with observations of Lya (see 
Lemoine et al. 1999 for a review). Linsky et al. (1993, 1995) found D/H = 1.60+o:i9 x !0~ 5 
toward Capella. Vidal-Madjar et al. (1998) reported evidence for a factor of ~ 2 difference in 
D/H between two components on the line of sight toward the white dwarf G191-B2B (d = 75 pc). 
Sahu et al. (1999) re-evaluated the GHRS data as well as new STIS echelle spectra of G191-B2B 
and concluded that, as a result of improved instrumental characterization (echelle scattered light 
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correction near Lya), D/H values in both components appeared to be consistent with the usual 
LISM value. Questions remain, however. For example, Howk & Sembach (2000) found a different 
STIS echelle scattered light correction, one that agrees with the GHRS Lya profile (Vidal-Madjar 
2000). The resolution of the conflicting conclusions about the G191-B2B sightline may be provided 
by observations of higher Lyman series lines (Ly/3 and Ly7) by the Far Ultraviolet Spectroscopic 
Explorer (FUSE, Moos et al. 2000). 

The Interstellar Medium Absorption Profile Spectrograph (IMAPS) provides a new window 
to study deuterium in the Galaxy by virtue of its wavelength coverage and very high spectral 
resolution (A/AA ~ 80,000) which alleviate many of the previous obstacles to accurate D/H 
measurements. Prior to the IMAPS deuterium studies, D/H measurements in our galaxy after 
the Copernicus mission have been confined to the local ISM, since the saturated core of H I Lya 
envelopes the D I Lya feature for 7V(H I) ^ 1 x 10 19 cm" 2 and higher Lyman lines lie beyond the 
reach of HST. The goal of the IMAPS deuterium program is to obtain high-quality measurements 
of D/H on sightlines toward bright OB stars beyond the LISM. In the first paper of this series, 
Jenkins et al. (1999a, hereafter Paper I) measured the D/H abundance ratio toward <50ri A using 
IMAPS spectra. They found D/H = 0.74^Q'{g x 10~ 5 , a value that is much lower than that found 
in the LISM. Paper I also showed that the low abundance of D toward 5 On A is not accompanied 
by an overabundance of N and O relative to H, as would be expected if the gas had been more 
thoroughly cycled through stellar interiors. We note that two other deuterium measurements using 
the Tubingen echelle spectrograph (resolution A/AA ~ 10,000) on the ORFEUS-SPAS //mission 
have been recently reported: Golz et al. (1998) found D/H = 0.8l°;4 x 10" 5 toward BD +28° 4211, 
an O-type subdwarf ~100 pc away, and Bluhm et al. (1999) measured D/H = l^t^ x 10~ 5 
toward BD +39° 3226, a sdO star located ~ 270 pc from the Sun. 

In this paper we analyze the lines of sight toward two well-known hot stars, 7 2 Velorum (HD 
68273) and £Puppis (HD 66811). They are among the brightest stars in the sky at 1000 A, well 
studied by Copernicus, have H I column densities <; 10 20 cm" 2 , and hence are prime D/H targets 
for IMAPS. These stars were observed under an ORFEUS-SPAS //Guest Investigator program. 
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Table 1. Stellar data 



Name 


HD 






Spectral 
Type 


V 

(mag) 


vsin i 
( km s _1 ) 


( km s _1 ) 


d h 

(PC) 


CPup 
7 2 Vel 


66811 
68273 


253? 98 
262? 81 


-4? 71 
-7? 70 


05 Iaf 
WC8+08 


2.25 
1.78 


211 


-18.04 
-17.48 


429+#° 
258lg 



a VLSR = Vhelio + Aw L5jR 

b Schaerer et al. (1997) 
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The sightlines to 7 2 Vel and £ Pup have been studied extensively (j 2 Vel: Morton & Bhavsar 
1979; Sahu 1992; Fitzpatrick & Spitzer 1994; (Pup: Morton & Dinerstein 1976; Morton 1978), 
including D/H measurements with observations from the Copernicus satellite. York & Rogerson 
(1976) found D/H= 2.0t 1 1 7 x 1(T 5 toward Vel. Vidal-Madjar et al. (1977) studied the sightline 
toward CP U P an d found that solutions consistent with the data permit a wide range of D/H 
values (1.7 x 10~ 5 ^ D/H ^ 2.5 x 10 -4 ). 7 2 Vel and CP U P are the most luminous stars in the 
Vela-Puppis region. Their role in powering the Gum Nebula and their relationship to the Vela 
OB2 association, Vela supernova remnants, and various structures in the ISM in Vela- Puppis 
have been analyzed by Sahu (1992). Table | summarizes some basic properties of the two stars. 
The distances listed were derived from HIP PARC OS data by Schaerer, Schmutz &: Grenon (1997), 
placing 7 2 Vel significantly closer than earlier estimates. 

7 2 Vel is the closest (and brightest) Wolf-Rayet star and consequently it has been studied 
intensively for more than a century (van der Hucht et al. 1981). Nevertheless, due to the 
complexity of the stellar system, its fundamental stellar parameters have been continuously 
debated and revised (as recently as 1997 when several papers were published, e.g., Schaerer et 
al. 1997; Schmutz et al. 1997). 7 2 Vel is a double- lined spectroscopic binary star composed of a 
WR star (WC8) and a late O star (08 III) of comparable brightness. The complex properties of 
this system are evident in the significant changes in the UV spectrum, and this behavior will be 
addressed in §^1] when we discuss our determination of the interstellar hydrogen column density 
through measurements of the Lya absorption feature. 

£Pup is an extremely luminous, massive star (04 Iaf) and the brightest O-type star in the 
sky. It is a prototype for O-star spectral properties and variability (Morton & Underhill 1977; 
Massa et al. 1995) and theoretical understanding of radiation-driven winds in massive stars (e.g. 
Puldrach et al. 1994). CP U P is believed to be a single star. 

The high velocity resolution provided by IMAPS opens the way for more accurate line profile 
measurements in the far ultraviolet spectra of bright stars. Our observations of 7 2 Vel and CP U P 
described in §[2] build on our earlier results for 5 Ori A described in Paper I. We follow with 
discussions of how we derived D I column densities (§||), H I column densities (§||), and values for 
D/H and N/H (§|5|). The paper concludes with a discussion in §|| of the possible implications of 
the differences in D/H for the three stars covered in this paper and Paper I. 

2. OBSERVATIONS 

The spectra analyzed in this paper were obtained with IMAPS during the ORFEUS-SPAS II 
mission (STS-80) in late 1996 (see Hurwitz et al. 1998). The design of IMAPS and its performance 
on earlier flights have been described in detail by Jenkins et al. (1988, 1996). In Paper I we 
described the in-flight performance of IMAPS during the 1996 mission. Here we only summarize 
the instrument's principal characteristics. 
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IMAPS is an objective-grating echelle spectrograph designed to record the far-UV spectrum 
of bright stars at high spectral resolution. The optical design consists of a wire grid collimator 
to reject off-axis light from stars other than the target, an echelle grating with a 63° blaze angle, 
and a parabolic cross disperser. The spectral format is imaged on a solid KBr photocathode, 
whose electrons are magnetically focussed on a windowless, back-illuminated CCD with a 320x256 
pixel format. The nominal wavelength range of 930-1150 A is obtained in four selectable tilts 
that span the free spectral range of the echelle grating. The gratings were coated with LiF over 
aluminum, providing excellent throughput longward of 1000 A. Although the reflectivity of LiF 
drops substantially shortward of 1000 A, IMAPS achieved a useful throughput even in the 930 
- 980 A region. The spectral resolution in IMAPS spectra obtained during the 1996 flight was 
approximately A/AA ~ 80, 000, or < 4 km s^ 1 . The telluric O I lines (e.g. A950.112 near Ly5 - see 
Fig H) have FWHM ~ 5 km s , but these lines are probably partly resolved in IMAPS spectra 
(see Jenkins et al. 1999b). 

7 2 Vel and £ Pup were chosen for this IMAPS Guest Investigator program because they 
were available in 1996 November, have a flux near Ly5 > 10~ 9 erg cm -2 s _1 A -1 , have 
vsini £ 100 km s -1 , and N(U I) £ 10 20 cm" 2 . 

7 2 Vel was observed by IMAPS on 1996 November 27 14:34-15:11 and 16:00-16:47 UT for a 
total of 4224 s. (Pup was observed on 1996 November 26 18:44-19:23 and 20:16-20:57 UT also for 
a total of 4224 s. The exposure time at each echelle grating position ranged from 817.6 s to 1226.4 
s. These exposure times were chosen to obtain spectra with good signal-to-noise (S/N) ratios near 
the cores of the D I Ly<5 (949.485 A) and Lye (937.548 A) lines. Typically, we found S/N=10-15 
in the local continuum near the Ly<5 and Lye interstellar features. 

3. COLUMN DENSITY OF ATOMIC DEUTERIUM 

3.1. General Considerations 

The data reduction and analysis of the 7 2 Vel and £ Pup spectra and the determination of 
the D I column densities were identical to those described in Paper I for the line of sight toward 
dOri A. IMAPS spectra were analyzed to determine the total iV(D I) (§ |3.3| and §p.4|). A large 
number of International Ultraviolet Explorer {IUE) high resolution Lya spectra were obtained 
from the National Space Science Data Center archive to determine the total H I column density 
(iV(H I) ) in §§. We note that given that N(R I) ~ 10 20 cm~ 2 on these sightlines and that D/H 
~ 10 -5 , the gas that gives rise to the H I Lya damping wings is the same gas responsible for the 
D I Ly5 and Lye features. 

In §[T^ we show that N I is a good tracer of H I and D I on these sight lines. The line of 
sight column density per unit velocity [N a (v), see Paper I] for N I, defined by a range of N I lines 
recorded in the IMAPS spectra, provided a velocity template for modelling the D I profiles. The 
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high S /N for the N I lines helped to constrain the model profiles that gave an acceptable fit to the 
lower S/N D I lines. This can prevent noise in the D I profiles from giving arbitrarily large or 
small N(D I) at specific velocities. We did not decompose the velocity profiles of N I or D I into 
separate (blended) Gaussian components - the D/H ratios determined here were based on total 
column densities for each sight line. 

Using the method outlined in Paper I, we corrected for the effects of the weak lines Fe II 
A937.652 and H 2 Lyman 14-0 P(2) A949.351. The Fe II line is located at +51 km s" 1 in -y 2 Vel 
and at +53 km s _1 in CP U P on the heliocentric velocity scale whose zero point is the laboratory 
wavelength of D I Lye (see Figs. [l| and |2|). The H2 line is located at +10 km s _1 in j 2 Vel and at 
+12 km s _1 in CP U P on the D I Ly<5 heliocentric velocity scale. The strength and shape of these 
features were computed from other transitions of the same species recorded at longer wavelengths 
(and higher S/N) in the IMAPS spectra of 7 2 Vel and CP U P- The computed central residual 
intensities of these features are nearly identical for the two stars: 0.66 for Fe II A937.6 and 0.95 
for H_2 A949.3. The other potentially contaminating lines noted in Table 1 of Paper I can safely be 
ignored. 

The possible contamination of the Ly<5 order by scattered light from an adjacent order that 
contains the strong absorption lines of the N I AA954 multiplet was evaluated in the same manner 
as described in Paper I. The x 2 analysis found that any residual contamination of the spectrum 
in the vicinity of Ly<5 by the pattern of saturated N I features was < 1% of the continuum in 
both 7 2 Vel and CP U P- Amplitudes for this correction larger than ~ 1% caused unacceptably bad 
deviations in the bottom of the H I Ly5 profile. Given the general noise characteristics of the Ly<5 
profiles, this is a negligible effect. 

For the case of 7 2 Vel, there is another potential contamination source. A companion star, 

7 1 Vel, is located 42" to the southwest of 7 2 Vel. Its light should be accepted along with that 
from 7 2 Vel, since IMAPS is an objective-grating spectrograph without an entrance slit to reject 
unwanted sources. (The grid collimator rejects light coming from more than 1° from the axis, 
however.) 7 1 Vel is 2.5 magnitudes fainter than 7 2 Vel in V and has a spectral classification 
of B1IV (Hofneit & Jaschek, 1982). Even though 7 1 Vel is cooler than 7 2 Vel, its flux is not 
diminished much below the peak of the Planck distribution at our wavelengths of interest, so we 
expect its intensity at a given wavelength to be not more than about a factor of 10 fainter than 

7 2 Vel. Fortunately, when the observations were taken the roll angle of the spacecraft about the 
optical axis of IMAPS (governed by the position of the Sun in the sky) was such that the faint 
spectrum of 7 1 Vel was displaced along the cross dispersion direction toward the long-wavelength 
part of the format. As a consequence, any spectral segment in the spectrum of 7 2 Vel had light 
superposed on it from shorter wavelengths in 7 1 Vel (the separation was slightly less than the 
distance between 3 echelle orders). The rapid decline in the sensitivity of IMAPS toward shorter 
wavelengths thus amplified the factor of 10 disparity in relative fluxes at any given position on the 
format. Therefore, the contamination of the 7 2 Vel spectrum by light from 7 1 Vel is negligible. As 
a final note, on visual inspection we are unable to see any ghost-like spectrum of 7 1 Vel on top of 
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the spectrum of j 2 Vel. 

The Lyman line profiles for j 2 Vel and £ Pup are shown in Figures |l] and ||. These profiles 
have been corrected for Fe II and H2 line absorption, as described above. The background levels 
near the D I lines were determined from the broad, saturated cores of the adjacent H I profiles. 
As explained in Paper I, we compared resolution-degraded forms of the IMAPS profiles with those 
recorded by Copernicus to test the proposition that the cores of the H I lines indeed represented 
the zero-intensity levels in the vicinity of the adjacent D I features. This was done to check that 
we were not being misled by an effect from possible strong, broad wings extending away from the 
main peak in the instrumental profile of IMAPS. We concluded that indeed the cores of the H I 
profiles provided very good estimates of the zero levels in the vicinity of the D I lines. 

The best answers for iV(D I) and the deviations permitted by the data were evaluated 
by minimizing x 2 ( see Lampton, Margon & Bowyer 1976 and Bevington & Robinson 1992 for 
details) when all of our unknown parameters were allowed to vary. The ten free parameters for 
7 2 Vel are (1) N(D I) , (2) the gas temperature T, (3-8) the continuum slopes, Y-intercepts, and 
background levels for Ly<5 and Lye, (9) a shift in the velocity zero point between the D I lines 
and the N I template, and (10) a coefficient for scaling the N I contamination signal in Ly<5 (see 
§3.2 in Paper I). The gas temperature is an important parameter because the D I line can be 
broadened significantly compared to the N I template by thermal motions; see §4.1 in Paper I for 
details. There are three additional free parameters for £ Pup because of the addition of Ly7 to 
the analysis. Note that the linear continuum fitting (with specific velocity limits as given in §§[3.3] 
and |3.4| ) is an integral part of the x 2 evaluation, i.e., the deviations of the continuum levels are 
considered, in addition to the behavior inside the D I profiles. The same holds for the zero level as 
defined by the bottom of the adjacent H I profile. We used Powell's method (Press et al. 1992, p. 
406) to find the minimum x 2 ■ We then set confidence limits by increasing (or decreasing) iV(D I) 
and T with the other parameters freely varying until x 2 increased by the appropriate amount for 
the confidence limit of interest with two useful parameters, N(D I) and T. We also used different 
initial values to establish that the minimum x 2 is unique. 
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Fig. 1. — Line profiles of LyS and Lye from IMAPS spectra for Vel. The zero-point of the velocity 
scale is computed with respect to the laboratory wavelengths of the D I lines. The narrow width 
of O I* A950.112 telluric line indicates the spectral resolution is A/AA ~ 80, 000. 
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Fig. 2. — Line profiles of Ly7, Ly<5, and Lye from IMAPS spectra for ("Pup. The zero-point of the 
velocity scale is computed with respect to the laboratory wavelengths of the D I lines. 
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3.2. Velocity Profile Templates 

In Paper I we discussed the benefits of obtaining a velocity profile based on high quality data 
for N I and O I, two species that are not significantly depleted in the ISM and which have very 
similar ionization balances to those of D I and H I (Ferlet 1981; York et al. 1983). This profile 
information is helpful in constraining the variety of possible interpretations that would produce 
acceptable fits with the D profiles. For N I, we used IMAPS spectra of the 10 lines in the multiplets 
at 952.4, 953.8, 954.1 and 1134.7 AQ. For all of the N I lines except those in the 952.4 A multiplet, 
the background level was easily defined because the stronger lines were saturated. Of course, these 
lines were useful only for defining the behavior of the gas at velocities somewhat removed from 
the line core. Since the lines in the multiplet at 952.4 A were not saturated, we had to determine 
the background level by a different method. For both 7 2 Vel and £Pup, there are Ul scans of this 
weak multiplet available in the archive of spectra recorded by the Copernicus satellite (Rogerson 
et al. 1973). After comparing our IMAPS spectrum degraded to the resolution of Copernicus with 
the actual Copernicus scans, we determined the adjustments to the background levels that were 
needed for the IMAPS spectra of this N I multiplet. The backgrounds caused by grating scatter in 
the Copernicus spectra were taken from the observed count rates in the bottoms of the hydrogen 
Ly7 and Ly5 lines, and these levels were subtracted off before the comparison was made. 

As was done for our analysis in Paper I, we adopted a method developed by Jenkins &; 
Peimbert (1997) to create a composite profile for the column density of N I as a function of 
velocity from the 10 lines in the four multiplets, with the weak lines in each case defining the main 
part of the profile and the strong lines outlining the exact shape of the profile's wings, well away 
from the saturated part of the line. For places where there was overlap in the useful portions of 
the lines, there was satisfactory agreement. The /-values of Goldbach et al. (1992) were adopted 
for our analysis. 

Figure ^ shows the N a (v) profiles derived for N I toward 7 2 Vel and ("Pup. In Paper I we 
showed that the profiles are unlikely to be contaminated by additional contributions from telluric 
absorption. iV(N I) is ~1% of N(0 I) in the Earth's atmosphere at the altitude of IMAPS (~305 
km) during the ORFEUS-SPAS II mission. Even for our strongest N I transition used to derive 
the far wings of N a (v) for N I (A1134.98), the telluric absorption should be only one third as 
strong as the O I* A950.112 feature present in Figures || and ^, and at a heliocentric velocity of 
+12.6 km s _1 it would be buried in the saturated portion of the profile. Consequently, telluric N I 
has a negligible effect on the interstellar N I profiles. 

O I is the best tracer of H I in the ISM since the ionization potential of O I (13.56 eV) is 
nearly identical to that of H I . The ionization potential of N I (14.53 eV) is only slightly greater 
that that of H I . Furthermore, both O and N are coupled to H by resonant charge exchange 



7 One of the lines in the 1134.7 A multiplet, the line at 1134.165 A, had to be omitted from consideration for £Pup 
because it was too close to the edge of the image format. 
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Fig. 3. — Profiles for N a (v) (thick solid lines) for N I toward j 2 Vel and £ Pup derived from IMAPS 
spectra. For 7 Vel, N a (v) for the wings of O I A1039.23, adjusted for the difference in cosmic 
abundance between O and N (—0.90 dex), is shown by the thin histogram line. The single data 
point and error bar correspond to N(0 I) from Fitzpatrick & Spitzer (1994) for O I A1355.6 (see 
text for details). The dashed lines show the expected shapes of the D I profiles for a thermal 
broadening for T = 6070 K (j 2 Vel) and T = 9550 K (£ Pup) of the respective N I profile favored 
by our most likely solutions in Tables S and y. 
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reactions. As discussed by Sofia & Jenkins (1998) and Jenkins et al. (2000), N I closely follows O I 
and H I unless n e 3> rani- 

In the IMAPS wavelength band there is no suitable set of O I transitions to completely define 
a velocity profile template. The available O I lines are highly saturated or are not detectable, 
either because they are buried in the core of a Lyman series line or because they are too weak. In 
Paper I we made use of an archival HST spectrum of the very weak intersystem transition of O I 
at 1355.6 A for 5 0ri A. There is no such spectrum available for £Pup, and the measurement of 
this line in the spectrum of 7 2 Vel taken by Fitzpatrick & Spitzer (1994) is too noisy to adequately 
define the shape of the O I profile. Therefore, we were forced to use the next best option, N I, to 
define the velocity profile template. 

Although we do not have a complete N a (v) profile for O I, we tested the assumption that 
N I traces H I by comparing selected portions of the O I N a (v) profile derived from available line 
profiles in 7 2 Vel. First, we computed N a (v) from the wings of O I A1039.230. Second, we used 
the GHRS spectrum of O I A1355.6 (Fitzpatrick & Spitzer 1994) to define N a (v) in the core of the 
line. The computed N a (v) was scaled by the difference in the N and O solar abundances (—0.90 
dex, Anders & Grevesse 1989). The results, shown in Fig ^, demonstrate that the wings and core 
of N a (v) for N I and O I are in excellent agreement. Consequently, we have high confidence that 
the N a (v) profile for N I provides an accurate model for the velocity distribution of H I . 



3.3. 7 2 Velorum 

The continuum near the Ly<5 line shown in Figure [l| was determined in the velocity range 
(—40, —10) km s _1 on the blue side of the D I line and (+170, +205) km s _1 on the red side. For 
Lye, the continuum limits were (—14, +8) km s _1 and (+170, +210) km s _1 . The background level 
was determined from the core of the two H I lines in the velocity range (+60, +140) km s _1 . The 
D I profiles were fit in the velocity range (—10, +39) km s _1 for Ly5 and (+8, +39) km s _1 for 
Lye. With a spacing of independent velocity samples of 1.25 km s _1 , this resulted in 276 degrees 



of freedom in the x analysis to fit the 10 free parameters described in §3.1. 



By adjusting our estimate for the noise, which was only known to limited accuracy beforehand, 
we achieved a minimum value for \ 2 equal to 247.0 at a column density of 1.12 x 10 15 cm -2 . 
For 276 degrees of freedom, there was a 90% chance that we would have found this value of x 2 
or greater. With this 90% confidence level, we arrived at a conservatively high estimate for the 
noise. It then followed that this noise level, which is perhaps higher than the real noise, gave 
a conservatively large confidence interval for N(D I) , as determined by how rapidly the values 
of x 2 deviated from the minimum value. With the noise level having been set in the manner 
just described, we explored for 90% and 99% confidence limits for N(D I) (i.e., 1.65a and 2.58a 
deviations), which correspond to x 2 ( m in)+4.6 and x 2 ( mm )+9-2, respectively. Table [2] lists iV(D I) 
and T for the best value and these limits. The ±90% confidence limits on the model D I profiles 
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(cross-hatched regions) for Ly<5 and Lye are compared with the observed deuterium profiles in 
Figure [|. 

In Figure || we also show with a dashed line the expected shape and depth of the D I profile for 
N(D I) =7.7 x 10 14 cm -2 . This corresponds to a D/H ratio of 1.5 x 1CT 5 , assuming the H I column 
density derived below in §||. For this case, x 2 ~ X 2 ( mm ) = 52.4, which is clearly an unacceptable 
fit. The dot-dashed line in Figure |] corresponds to the D I profile for N(D I) =3.8 x 10 14 cm~ 2 , 
the D I column density necessary to achieve D/H=0.74 x 10 -5 , the value found toward 6 Ori A in 
Paper I. A similar demonstration for the measurement of iV(H I) is be given in §|5| 



3.4. CPuppis 

The measurement of N(D I) toward CP U P used the same methodology described in the 
previous section with the following modifications. In addition to Ly<5 and Lye we were able 
to include the D I profile of Ly7 in the x 2 analysis. The blue and red regions for the linear 
continuum fits were (—40,-10) km s _1 and (+170, +190) km s _1 for Ly7, (—39,-10) km s _1 
and (+157, +190) km s" 1 for Ly<5, and (-4, +3) km s" 1 and (+180, +220) km s" 1 for Lye. The 
background levels were determined from the saturated cores of the H I lines over these velocity 
limits: (+55, +130) for Ly7 and Ly<5, (+65, +115) for Lye. The model D I profiles were fit over 
the range (—5, +40) km s _1 for Ly7, (0, +40) km s _1 for Ly<5, and (+3, +37) km s _1 for Lye. 
This resulted in 426 degrees of freedom to simultaneously fit the 13 free parameters described in 
§[3+| The minimum x 2 ls 388.5 and corresponds to the most probable value of N(D I) for £Pup, 
namely 1.30 x 10 15 cm" 2 . The results of the x 2 analysis are summarized in Table |||, including the 
90% and 99% confidence limits. The same conservative treatment of the noise described for 7 2 Vel 
was followed for £ Pup. Figure |H| shows the observed D I Lyman line profiles and the range in the 
best fit model profiles and continua fits allowed by the 90% confidence limits. 
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Table 2. Limits for 7V(D I) Toward 7 2 Vel 





N(B I) 


T 




Significance 


(10 15 cm' 2 ) 


(K) 


x 2 


Minimum N(D I) at the 99% confidence limit 


0.96 


6630 


255.6 


Minimum N(D I) at the 90% confidence limit 


1.00 


6530 


251.6 


Best N(D I) 


1.12 


6070 


247.0 


Maximum iV(D I) at the 90% confidence limit 


1.27 


5280 


251.7 


Maximum iV(D I) at the 99% confidence limit 


1.34 


4870 


256.1 



Table 3. Limits for JV(D I) Toward C Pup 





N(B I) 


T 




Significance 


(10 15 cm" 2 ) 


(K) 


x 2 


Minimum iV(D I) at the 99% confidence limit 


1.07 


10600 


397.7 


Minimum iV(D I) at the 90% confidence limit 


1.13 


10400 


393.1 


Best N(B I) 


1.30 


9550 


388.5 


Maximum N(D I) at the 90% confidence limit 


1.49 


8500 


393.1 


Maximum iV(D I) at the 99% confidence limit 


1.58 


7900 


397.7 
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Fig. 4. — Observed and model line profiles of D I Ly5 and Lye for 7 2 Vel showing the region 
bounded by the 90% confidence limits on iV(D I) (cross- hatched region), derived from both lines 
collectively. Assuming the value of iV(H I) derived for 7 2 Vel in §||, the dashed line illustrates the 
D I line shapes necessary to achieve D/H = 1.5 x 1CP 5 and the dot-dashed line corresponds to the 
D I profiles necessary to obtain D/H = 0.74 x 10~ 5 (5 Ori A, Paper I). 




Fig. 5. — Observed and model line profiles of D I Ly7, Ly<5, and Lye for CP U P showing the region 
bounded by the 90% confidence limits on iV(D I) (cross-hatched region), derived from all three 
lines collectively. Note that the preferred continuum levels near Lye change for different values of 
N(D I) . 
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4. COLUMN DENSITY OF ATOMIC HYDROGEN 

4.1. Analysis Methodology 

The primary objective of this IMAPS program is to determine D/H ratios of sufficient 
accuracy to test for spatial inhomogeneities. To achieve this goal, we must strive for a precision 
in N(H I) that is as good as (or preferably better than) that of the deuterium measurement. We 
must also understand the magnitude and nature of the uncertainties in iV(H I) and iV(D I) so that 
we can realistically assess the uncertainty in D/H for comparison to other measurements in the 
ISM. For these reasons, we decided to conduct our own investigation of the H I column densities 
toward 7 2 Vel and £Pup, even though a number of measurements of iV(H I) toward these stars 
are already available in the literature (e.g. Jenkins 1971; York & Rogerson 1976; Vidal-Madjar et 
al. 1977; Bohlin, Savage, & Drake 1978; Shull k Van Steenberg 1985; Diplas & Savage 1994). 

We described our method for determining iV(H I) in Paper I; here we provide a summary. 
In §4.2 and §4.3 we discuss details particular to the CP U P an d 7 2 Vel sight lines and present the 



results. We used the Lya line to determine the total N(R I) ^\ Since H I Ly7 and Ly<5 are on the 
flat part of the curve of growth, high velocity gas could influence their profiles. Lya is immune 
from this potential problem. Due to the great breadth of the Lorentzian wings (±1000 km s _1 
- see the velocity scale in Figure p|), Lya absorption due to any high velocity interstellar gas 
(±100 km s" 1 ) that may be present is confined to the saturated core of the profile. In principle, 
some high velocity H I , which is too widely dispersed in velocity to show up in the D I profiles, 
could contribute to the Lya damping wings. However, this gas is not detected in the strong lines 
N I A1134.98 and O I A1039.23 observed with IMAPS. It is clear that such high velocity gas, 
if present, would have a low column density and a negligible effect on the Lya damping wings. 
Since the Lya line is not covered by IMAPS, was not observed by the HST spectrographs, and 
the available Copernicus Lya scans suffer from several problems (see Paper I), we analyzed high 
dispersion IUE observations of Lya to determine iV(H I) . Furthermore, these stars were observed 
many times with IUE over the course of many years, and this offers an opportunity to evaluate 
potential sources of systematic error. For example, 7 2 Vel is a spectroscopic binary, and the large 
IUE database allowed a search for systematic changes in the derived N(H I) as a function of 
orbital phase. 

After screening and retrieving the IUE data^ of interest and correcting for interorder scattered 
light (see Paper I), we used an approach first used by Jenkins (1971) to constrain the H I column 



8 The Copernicus studies of D/H on the 7 2 Vel and CP U P sight lines used the much weaker damping wings of Ly/3 
or Ly7 to determine JV(H I) . 

9 We used the standard IUESIPS calibration of the IUE high dispersion spectra instead of the NEWSIPS calibration 
to avoid potential problems with background corrections and zero-level determination near Lya found in some IUE 
SWP high dispersion spectra processed with NEWSIPS (see Massa et al. 1998), which could adversely affect our 
measurements. 
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density, i.e., we determined the iV(H I) which provides the best fit to the Lya profile with the 
optical depth r at a given wavelength A calculated from the expression 

r(A) = JV(H I)<r(A) = 4.26 x 10~ 20 N(H I)(A - A )~ 2 (1) 

where Ao is the Lya line center at the velocity centroid of the hydrogen (note that the useful 
portion of the Lya profile is entirely due to the Lorentzian wings and the effects of instrumental 
and Doppler broadening can be neglected). As was done for D I , we also determined the important 
free parameters which can be adjusted to fit the H I Lya absorption profile, then we found the 
set of parameters that minimized x 2 using Powell's method. We set confidence limits on the H I 
column as described in §3.1 with only one parameter of interest, iV(H I) . The other (uninteresting) 
free parameters we selected for fitting the H I profile were three coefficients which specify a 
second-order polynomial fit to the continuum, and an additive correction to the flux zero point. 
The continuum was fit to the spectrum in several windows covering the range 1185 A to 1276 A. 
Despite our use of the Bianchi & Bohlin (1984) correction for interorder scattered light, in many 
cases inspection of the flat-bottomed, saturated portion of the Lya profile showed that the zero 
intensity level was not quite correct. A zero point shift in the flux scale, one of the terms for 
evaluating x 2 > was determined from a region within the saturated core. For both stars, the zero 
point of the Lya wavelength scale was set by placing the N I A 1200 multiplet in agreement with 
the IMAPS N I profiles. The Lya profile was then fit only to the red wing because of the presence 
of strong stellar features superposed on the blue wing. However, the uncontaminated portion of 
the blue wing of Lya was checked for consistency with the fit to the red wing and found to be in 
good agreement. 

4.2. £ Puppis 

In Paper I, considerable attention was paid to the spectroscopic binary nature of 5 Ori A 
and whether or not this could cause systematic errors in A r (H I) . A similar analysis of 7 2 Vel is 
presented in § |4.3| . The determination of N(H I) toward CP U P was l ess complex than the other 
stars observed by IMAPS to study D/H. As far as is known, CP U P is n °t a spectroscopic binary 
star. It is an 04 supergiant, so the stellar H I Lya absorption line makes a negligible contribution 
to the H I absorption profile. The star is possibly a non-radial pulsator with very weakly variable 
stellar absorption lines with a period of 8.5 hr (Reid & Howarth 1996), and the variability of 
the stellar wind P-Cygni profiles is well known and studied (e.g., Prinja et al. 1992; Howarth, 
Prinja, & Massa 1995) and shows significant power at 19.2 hr and 5.2 day periods. Both of these 
sources of spectral variations are expected to have little impact on the interstellar H I column 
density derived from Lya, but this can be checked given the large number of observations and 
good temporal sampling. 

There are more than 200 high-dispersion observations of CPup i n the IUE archive, primarily 
because the star was selected for intensive IUE observing programs to study stellar wind variability 
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in massive stars. We concentrated on the spectra obtained in three specific multi-day observing 
runs in 1986, 1989, and 1995. In 1995 CP U P was observed continuously for 16 days (dubbed the 
"MEGA" campaign, Massa et al. 1995). We omitted three observations from these observing 
sessions because the archival data were corrupted or unavailable. This left 189 observations for 
measurement of iV(H I) . All of the observations were obtained with the IUE large aperture and 
the signal-to-noise ratios are comparable; the CP U P data set is more uniform than the IUE Lya 
data used for measuring iV(H I) toward 5 Ori A and 7 2 Vel. 

The iV(H I) values derived from each individual observation are plotted in Figure ||, along 
with the la uncertainties, versus the observation date. Table |] summarizes the mean H I column 
density <iV(H I) > and the root mean square (rms) dispersion a obtained from the 1986, 1989, 
and 1995 data sets analyzed separately, and <iV(H I) > for each data set is indicated with a heavy 
dashed line in Figure ||. We also list in Table || the formal error in the mean e = a/\f~N (here N is 
the number of measurements) and the reduced x 2 f° r the three data sets. 

In Paper I we discussed the potential sources of systematic error in the derivation of iV(H I) 
from IUE spectra, and we suggested that the real uncertainty in the mean is likely to be greater 
than the formal estimate. The very large Lya data set for CP U P shown in Figure ^ appears to 
confirm that this is indeed the case for this star. From this figure one can see that the column 
densities derived from the 1995 data are systematically lower than the column densities derived 
from the 1989 observing run. The mean of the 1995 data set is lower than the mean of the 1989 
data by 0.54 x 10 19 cm -2 (0.0256 dex), and this difference is substantially greater than the formal 
error estimates (e). 

This systematic error could be the result of secular changes in the stellar mass loss or 
circumstellar environment, a previously unrecognized stellar variability with a long period, or 
instrumental effects. Or perhaps some aspect of the interstellar sight line changed during this 
period. Since the IMAPS observations for N(D I) were made almost two years after the MEGA 
campaign, it is possible that a similar systematic error is present in any IUE estimate of iV(H I) 
that we choose for comparison with the IMAPS N(D I) . The exact magnitude of this systematic 
uncertainty is difficult to estimate, but we find that iV(H I) obtained from additional large 
aperture spectra of £ Pup taken at other epochs are in good agreement with the estimates in Table 
||. This indicates that the systematic error is not much larger than 0.0256 dex. For example, 
large aperture observations taken in 1988 May, 1989 December, and 1991 March, combined 
with the three <A r (H I) > values for the three data sets in Table ||, yield a mean N(H I) = 
9.34 ± 0.67 x 10 19 cm~ 2 . We conclude that 0.0256 dex is a conservative estimate of the la 
uncertainty in JV(H I) . We therefore adopt the unweighted mean of the three <iV(H I) > values 
in Table || as the best estimate of A r (H I) toward CP U P an d 0.0256 dex as its la uncertainty: 
N(B I) = 9.18 ± 0.54 x 10 19 cm" 2 . We elected to use a mean of the <iV(H I) > values in Table § 
instead of a straight mean of all of the individual measurements shown in Fig. || because the latter 
would give excessive weight to the 1995 data set due to the much larger number of measurements 
obtained during that observing campaign. 
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Fig. 6. — H I column densities toward C? U P derived from the 189 high-dispersion IUE observations 
obtained in 1986, 1989, and 1995, plotted versus the observation date. Each measurement is shown 
with ±lc7 error bars. All observations were obtained with the large aperture and have comparable 
signal-to-noise ratios. Each data set is separated by a vertical dotted line (i.e., there is a break 
in the x-axis at each vertical dotted line), and the mean H I column densities derived from the 
three data sets are indicated with heavy dashed lines. There appears to be a systematic difference 
between the column densities derived from the 1995 data and the earlier observations. 
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Table 4. Interstellar JV(H I) Toward CPup 



Data 


No. 


< iV(H I)> 




e c 


2 d 


Set a 


Spectra 


( cm -2 ) 


( cm -2 ) 


( cm -2 ) 




1986 


11 


9.23xl0 19 


0.47xl0 19 


0.14xl0 19 


0.69 


1989 


31 


9.42xl0 19 


0.43xl0 19 


0.08xl0 19 


0.56 


1995 e 


147 


8.88xl0 19 


0.46xl0 19 


0.04xl0 19 


0.69 


Adopted value f 


189 


9.18xl0 19 









a All A^(H I) measurements were derived from IUE high dispersion 
spectra using the large aperture. We have used only the observations 
obtained during three observing runs in 1986, 1989, and 1995. This table 
shows the mean H I column density, <A r (H I) >, derived from the 1986 
data only, the 1989 data only, and the 1995 data only. 

h a = rms dispersion (both a and <A r (H I) > were weighted inversely 
by the variances of the individual A^(H I) measurements). 

c e = error in the mean = c/(No. measurements) 0,5 . 

d Reduced % 2 = x 2 /(degrees of freedom), where 
X 2 = Ei { Wi (H I)- < 7V(H I) >] jo [AT(H I)] J 2 . 

c The IUE MEGA campaign (Massa et al. 1995). 

f Final adopted value of 7V(H I) is the average of <A r (H I) > for the 
three data sets. See text for details. 



-23- 



4.3. 7 2 Velorum 

As noted in §1, 7 2 Vel is a complex stellar system with a strongly variable UV stellar spectrum. 
Figure ^ shows two examples of the high-dispersion IUE spectra of 7 2 Vel that we have used to 
derive JV(H I) . These examples were selected to illustrate the variability of the P-Cygni emission 
lines in the vicinity of Lya, which are weakest at orbital phase ~0.5. Figure |7| also shows the 
complex and variable stellar absorption superposed on much of the blue wing of the interstellar 
Lya profile. This variability was a source of concern for this paper. Can we derive reliable H I 
column densities despite the complex spectral changes occurring in this star? Due to the stellar 
absorption, the blue wing of Lya was not useful for constraining iV(H I) . Fortunately, the red 
wing of Lya was clean and stable as a function of orbital phase. This can be seen by comparing 
panel (a) to panel (b) in Figure [?]. In this paper, we used only the red wing for fitting the Lya 
profile, and the resultant fits generally look quite good (see Figures [7| and ^). 

We examined iV(H I) as a function of orbital phase of the spectroscopic binary to check for 
systematic errors in the derived interstellar H I column density due to the complicated variability 
of the stellar spectrum. There are considerable discrepancies in the orbital elements published by 
various groups (Stickland & Lloyd 1990; Schmutz et al. 1997, and references therein). We adopt 
the orbital elements derived by Schmutz et al. (1997) who find a period of 78"?53 ± 0^01 d with 
velocity semiamplitudes of -f^wR = 122 ± 2 km s _1 and Kq = 38.4 ± 2 km s _1 . In addition to the 
usual intrinsic variability of WR stars, it is believed that the 7 2 Vel spectrum may also change as a 
result of periodic absorption of the O star component by the WR wind (Stickland & Lloyd 1990). 

After screening the IUE high-dispersion observations of 7 2 Vel, we were left with 42 spectra, 
34 small aperture observations obtained early in the IUE mission and 8 later observations obtained 
with the large aperture. In Figure |8] we plot the H I column densities derived from all of the 
observations as a function of the phase of the spectroscopic binary, and Table || summarizes the 
<N(H I) >, a, and e derived from the large aperture data only, the small aperture data only, and 
all of the data combined. While one might argue that a trend is apparent in Figure |§] with a 
minimum in the derived iV(H I) values at orbital phase ~ 0.5, this is a marginal result at best, 
and we do not believe that it should be taken too seriously. On the contrary, the generally good 
agreement of the H I column densities at various phases despite the dramatic changes in the stellar 
spectrum (see Figure 0) is encouraging, and we conclude from Figure || that the IUE data provide 
a good determination of <iV(H I) >. Since the small and large aperture measurements are in 
good agreement, we adopted the results from the combined data set. Systematic errors probably 
cause the true uncertainty in <iV(H I) > to exceed the formal error estimate provided in Table ||. 
Drawing from our experience with the more extensive data set for £Pup, we conservatively adopt 
0.0256 dex as the la uncertainty in iV(H I) toward 7 2 Vel as well: N(H I) = 5. 13 ±0.30 x 10 19 cm" 2 . 
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Fig. 7. — Samples of the IUE observations of 7 2 Vel used to derive the interstellar N(H I) : 
(a) SWP6175, spectroscopic binary phase = 0.33, and (b) SWP4719, phase = 0.48 (phases were 
calculated using the period and To from Schmutz et al. 1997). These spectra have been smoothed 
with a 5-pixel boxcar for display purposes only; the unsmoothed data were used to constrain N(H I) 
as described in the text. The best-fitting H I profiles (dotted lines) and continua (dashed lines) 
are overplotted on the data. Note the dramatic variability of the P-Cygni profiles. Note also 
the presence of complicated absorption structure which makes the blue wing of the Lya profile 
difficult to use for constraining iV(H I) . Both of these observations were obtained with the small 
IUE aperture. 
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Fig. 8. — Interstellar H I column densities derived from the 42 IUE observations of 7 2 Vel, plotted 
versus phase of the spectroscopic binary. Each column density is plotted with ±lcr error bars. Data 
obtained with the large aperture are indicated with open circles while small aperture results are 
shown with filled circles. 
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Table 5. Interstellar iV(H I) Toward 7 2 Vel 



Data 


No. 


<N(R I)> 




e c 


2 d 
Ay 


Set a 


Spectra 


( cm -2 ) 


( cm -2 ) 


( cm -2 ) 




Lg. Aperture 


8 


5.35xl0 19 


0.51xl0 19 


0.18xl0 19 


1.15 


Sm. Aperture 


34 


5.10xl0 19 


0.48xl0 19 


0.08xl0 19 


1.71 


All Data 


42 


5.13xl0 19 


0.49xl0 19 


0.08xl0 19 


1.62 



a All N(H I) measurements were derived from IUE high dispersion 
spectra. This table shows the mean H I column density, <A r (H I) >, 
derived from the large IUE aperture data only, the small aperture data 
only, and all of the data combined. 

b a = rms dispersion (both a and <iV(H I) > were weighted inversely 
by the variances of the individual iV(H I) measurements). 

c e = error in the mean = a /(No. measurements) 0,5 . 

d Reduced = X 2 / (degrees of freedom) , where \ 2 = 

Ei { [Ni (H I)- < N(R I) >] /a [N(H I)].} 2 . 
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5. DEUTERIUM AND NITROGEN ABUNDANCE RATIOS 

5.1. D/H 

Combining the results for D I and H I column densities derived in the previous sections, we 
determined the atomic D/H abundance ratio for 7 2 Vel and £Pup. Toward 7 2 Vel we found that 
D/H = 2.18j^Q32 x (the errors are 90% confidence limits). This result is a large departure 
from the value of D/H usually attributed as "typical" in the ISM of the Galaxy (~ 1.5 x 1CP 5 ). 
A large value for the average D/H ratio implies that even larger values could exist in individual 
components if other components have lower D/H values closer to the LISM ratio. Although the 
previous D/H measurement for 7 2 Vel (= 2.0^'^ x 10~ 5 , York & Rogerson 1976) appears to be 
in close agreement with the IMAPS result, we consider this to be coincidental in view of the 
magnitude of the uncertainties in the Copernicus measurement. 

For CPup we find JV(D I) = 1.30lg;$ x 10 15 cm 2 , a value slightly below the lower limit 
derived from Copernicus spectra by Vidal-Madjar et al. (1977). They found that a large range in 
N(D I) (1.5 - 20. x 10 15 cm" 2 ) was possible for this sight line due to its complexity and a lack of 
adequate constraints. Since the value of N(R I) derived by Vidal-Madjar et al. (1977) for C? U P is 
consistent with ours, their large range in D/H must be a result of the JV(D I) uncertainty. Much 
more is known now about the complexity of this sightline than at the time of the Copernicus 
study (e.g. Welty, Morton & Hobbs 1996). More importantly, the IMAPS spectra have sufficient 
spectral resolution that the structure of the £ Pup sightline in neutral hydrogen can be defined by 
the N a (v) profile for N I and used to constrain the determination of iV(D I) . Thus with higher 
spectral resolution and better knowledge of the velocity structure along the sightline, a more 
accurate value of iV(D I) and D/H is derived. We find D/H = 1.42±g;|| x 10~ 5 . The JV(D I) , 
N(R I) , and D/H results for the three stars studied by IMAPS are summarized in Table |(| where 
all errors are presented as 90% confidence limits. 

We now address a critical question, namely, is there sufficient uncertainty in the 7V(H I) 
and iV(D I) measurements to reconcile the 7 2 Vel D/H abundance ratio with the general result 
observed in the local ISM or with the D/H derived in Paper I for the sight line to 5 Ori A? Figure 
||a shows with dotted and dashed lines the expected appearance the Lya profile would have if 
our measurement of N(D I) is correct but the iV(H I) were made high enough so that D/H 
= 1.5 x 10" 5 . In Fi gure ^, the dotted line represents the H I profile that would be needed for 
the most probable value for N(D I) , while the dashed lines represent the 90% confidence limits. 
For comparison, we also show in Figure [9|a the best-fitting H I profile (solid line) derived in the 
previous section for this particular observation. Figure |||b is an analogous plot with N(R I) forced 
to take on a value so that D/H = 0.74 x 10 -5 as observed toward 5 Ori A (Paper I). From Figure 
H it is clear that forcing N(H I) to make D/H = 1.5 x 10~ 5 , assuming A^(D I) from § |3.3| , produces 
a very poor (and unacceptable) fit to the H I profile, and the fit becomes ridiculous if D/H 
= 0.74 x 10~ 5 . We have assumed the best-fit continuum (solid line) to compute the Lya profiles 
shown in Figure |]. We have also explored whether or not the D/H ratios can be reconciled by 
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Tablc 6. Abundance Ratios' 



7 2 Vel (Pup ,5 Ori A 



JV(DI) (10 15 


cm' 2 ) 


i 19+0.15 

i - iz -0.12 


1 so +0 - 19 

1 -* JU -0.17 


1 lfi+ - 29 

1 - 1D -0.20 


N(N I) (10 15 


cm -2 ) 


4.10 ±0.34 


7.62 ± 0.84 


6.19 ±0.51 


iV(Hl) (10 19 


cm -2 ) 


5.13 ±0.50 


9.18 ±0.89 


15.6 ± 1.52 b 


D/H (10- 


- 5 ) 


9 1O+0.36 


1 42+ - 25 


n 7/1 +0.20b 
u - ' 4 -0.15 


N/H (10" 




7.99 ± 1.02 


8.30 ± 1.22 


3.97 ±0.51 


D/N 




0.27 ±0.04 


0.17 ±0.03 


0.19 ±0.05 



a All errors presented in this table are 90% confidence limits 
(1.65a). 

b Errors in N(H I) for 5 Ori A from Paper I were revised 
to reflect the fact that potential systematic effects in N{H I) 
identified in this paper for CP U P rnay dominate the statistical 
errors quoted in the earlier paper. The la error in iV(H I) was 
set to 0.0256 dex. The errors in D/H for 5 Ori A are essentially 
the same as found in Paper I since they are dominated by the 
errors in iV(D I) . 
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Fig. 9. — Expanded plots of the H I Lya profile in 7 Vel shown in Figure 0b (SWP4719), 
again smoothed with a 5-pixel boxcar for display purposes only. The best-fitting profile and its 
corresponding continuum are indicated with solid lines. The broader damping profiles outside the 
best fit show the appearance of the profile when iV(H I) is forced to take on a value so that (a) 
D/H = 1.5 x 10~ 5 as generally observed in the local ISM, given the values of N(D I) in Table || 
[the dotted line assumes the most probable N(D I) while the dashed lines use the 90% confidence 
limits on N(D I)], and (b) D/H = 0.74 x 10~ 5 as observed toward 5 Ori A in Paper I. Note that 
the region between 1207 and 1210 A is sometimes affected by P-Cygni emission (which happens to 
be weak in this observation, see Figure 0). 
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adjusting the continuum placement, and we found that absurd continuum placements and shapes 
must be assumed to make D/H = 1.5 x 10~ 5 or (even worse) 0.74 x 10 -5 . 

As demonstrated at the end of §3.3, the changes in N(D I) needed to make D/H = 1.5 x 10 -5 
or = 0.74 x 10 -5 , assuming our most probable value of N(H. I) , gave a large mismatch with the 
data, which were clearly unacceptable. We conclude that spatial inhomogeneities in D/H in the 
ISM within ~500 pc of the Sun have high statistical significance, especially when contrasting the 
D/H ratios in the directions of 7 2 Vel and <50ri A, which were measured with the same instrument 
and techniques. 



5.2. Nitrogen Abundances 

The N I column density N(N I) was computed by integrating the column density profiles 
shown in Figure ^. The N(N I) results for the three targets are listed in Table ||. This table also 
contains the resulting N/H and D/N abundance ratios. We examined the potential error in iV(N I) 
by examining the errors in the portions of the various N I profiles used to construct N a (v). The 
la uncertainty in iV(N I) for -y 2 Vel and 5 Ori A is conservatively estimated to be 5% since the N I 
profiles were all of high quality. There may be some question about the structure in the core of 
N a (v) for £Pup. The uncertainty in iV(N I) for CP U P was estimated by supposing the spike in 
N a (v) at v ~ 19 km s _1 is due largely to noise fluctuations in the core of N I A952.5, the weakest 
N I line we detect on this sightline. If we truncate the N a (y) profile at the level of the secondary 
maximum at v ~ 24 km s^ 1 \N a [v)= 5.6 x 10 14 cm _2 ( km s -1 ) -1 ], we find that the area of the 
spike is 3.6 x 10 14 cm -2 , or 4.7% of the total. We then suppose that the total uncertainty in 
iV(N I) for CP U P is times this, or la = 5.1 x 10 14 cm -2 , to allow for the possibility that there 
may be other such fluctuations. 



6. DISCUSSION 

There are two principal conclusions of the IMAPS D/H program. First, the atomic D/H ratio 
in the ISM, averaged over path lengths of 250 to 500 pc, exhibits significant spatial variability. 
Differences in the atomic D/H ratio on long pathlengths in the ISM have been suspected for 
many years (Vidal-Madjar et al. 1978; Vidal-Madjar &: Gry 1984), but not until now have data of 
sufficient quality been available to evaluate and reduce statistical and systematic errors to levels 
where these differences are unequivocal. Second, we find no support for the simple picture that 
variations in D/H anticorrelate with those of N/H, i.e., one measure of how much the gas has 
been processed through stellar interiors. Figure |l0| shows the relationship between D/H and N/H 
for the three sightlines studied by IMAPS plus the white dwarf G191— B2B (Vidal-Madjar et al. 
1998; Sahu et al. 1999). We point out that some elements are systematically removed from the gas 
phase as they are incorporated into interstellar dust (Savage k, Sembach 1996), but the abundance 
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Fig. 10. — N/H abundance ratio as a function of D/H for the three stars studied by IMAPS plus 
the white dwarf G191-B2B. The abundance ratios for G191-B2B represent the full range of values 
reported by Vidal-Madjar et al. (1998) and Sahu et al. (1999). The average value of N/H in the 
ISM found by Meyer, Cardelli, & Sofia (1997) is indicated by a dashed line. The errors bars are 
90% confidence limits (1.65tr). 
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of N does not seem to be appreciably altered by this effect (Meyer et al. 1997). 

Beyond the effects from depletions onto dust, spatial variations in interstellar gas abundances 
can arise as a natural consequence of galactic chemical evolution and the changing influences of 
different stellar populations. The lack of an anticorrelation between N/H and D/H depicted in 



Fig. 1C indicates that the variability of D/H is not just a consequence of different mixing ratios 
of material with differing levels of stellar processing, as we might anticipate, for instance, from 
the variable addition of metal-poor, infalling gas from the Galactic halo (Meyer et al. 1994). 
We may need to go further and draw a distinction between contributions from stars that simply 
destroy deuterium and those that both destroy deuterium and enrich the medium with additional 
nitrogen. That is, we could envision some stars cycling material only through their shallow layers 
that are only hot enough to burn deuterium, while others eject material from much deeper layers 
where the synthesis of heavier elements has taken place. This additional level of complexity could 
explain the behavior that we observed. 

Global models of Galactic chemical evolution (Audouze & Tinsley 1974; Tosi 1988a, b; 
Dearborn, Steigman & Tosi 1996; Scully et al. 1996; Tosi et al. 1998) describe the destruction of 
D during stellar formation, evolution, and eventual mass loss. These models predict variations 
in D/H, N, and O abundances that are manifested as abundance gradients on a scale of ^1 
kpc. However, the predicted trends in galactic abundances may not accurately represent what 
is observable in the diffuse ISM. Tenorio-Tagle (1996) showed that the chemical enrichment of 
the diffuse ISM by OB associations, including supernovae from massive stars, is a slow process. 
Following a supernova explosion, chemically enriched ejecta remain clumpy and not well mixed 
with the diffuse ISM it encounters until it is incorporated into new star forming regions. This is 
a result of very long time scales for diffusion between different parcels of gas in the warm (10 4 K) 
and cold (10 2 K) phases of the ISM. 

The diffusion time scale for enriched gas to thoroughly mix with the warm diffuse ISM 
can be very long (tp > 10 10 yr - Tenorio-Tagle 1996). On the spatial scale sampled by the 
IMAPS observations, different sight lines may encounter regions with very different dynamical and 
chemical histories. Differential galactic rotation and random cloud motions are expected to stir the 
diffuse ISM and chemically enriched parcels of gas, but these parcels retain their distinct chemical 
properties until they are disrupted through photoevaporation, most likely by the formation of 
new massive stars. Tenorio-Tagle finds that diffusion is efficient only for the hot phase of the 
ISM (tf) < 10 6 yr), which accounts for a very small fraction of the total diffuse ISM. Only after 
the enriched gas and diffuse gas are highly ionized would the chemically enriched gas from the 
earlier generation of stars quickly diffuse into the ambient ISM. Thus, the time scale for mixing 
interstellar gases with different processing histories can be much longer than the chemical evolution 
time scale. However, it is possible that interstellar turbulence and its secondary phenomena 
may accelerate the mixing rate. Since the distribution of star forming regions (OB associations) 
shows large inhomogeneities on scales <^ 1 kpc, their corresponding chemical enrichment of the 
ISM may be expected to be nonuniform as well. This is perhaps revealed indirectly by variations 
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in H II region abundances (Peimbert 1999) and solar-type stars at similar galactocentric radii 
(Edvardsson, et al. 1993). 

Several processes unrelated to stellar nucleosynthesis may, under the right circumstances, 
alter the atomic D/H ratio of some parcels of interstellar gas (see Lemoine et al. 1999 for a 
review). D may be incorporated into HD (Watson 1973), but the fraction of molecular gas on our 
sight lines is very low (see Paper I). Differential radiation pressure on D and H (Vidal-Madjar 
et al. 1978; Bruston et al. 1981) may lead to a separation of D in some clouds near strong 
radiation fields. Adsorption of D onto dust grains (Jura 1982) may deplete D from the gas phase. 
Bauschlicher (1998) found that reactions of H and D with polycyclic aromatic hydrocarbon (PAH) 
cations might systematically provide some D enrichment in PAHs. A very different perspective has 
been offered by Mullan & Linsky (1999), who suggested that significant quantities of D may be 
formed in stellar flares from M dwarf stars and ejected into interstellar space. Some or all of these 
processes could be at work in the diffuse ISM and alter the atomic D/H ratio on individual sight 
lines, independent of the degree of chemical enrichment from stellar evolution. However, they have 
not yet been demonstrated to be quantitatively significant. Observational and theoretical tests of 
the efficiency and applicability of these processes are needed to better understand the mechanisms 
affecting the D/H ratio in the diffuse ISM. 

While the D/H ratios derived from IMAPS spectra are in the general range expected from 
galactic chemical evolution models, a factor of three variation in the mean D/H ratio on path 
lengths of several hundred pc is unexpected. The apparent lack of an anti-correlation of the D/H 
abundance ratio with the metallicity of the gas and its variability on smaller than expected scales 
suggests that other processes in the Galaxy may be masking more general chemical evolution 
trends. This may pose a problem for deriving a "primordial" D/H by extrapolating back from 
D/H measurements in the Milky Way to extragalactic absorbers at higher and higher redshifts, 
or to even to evaluate "primordial" D/H directly from high redshift observations, until we 
understand the reasons for these differences. The spatial variations found in this study underline 
the importance of high-quality D/H determinations. D/H measurements in more distant regions 
of the Galaxy are needed to determine whether the properties of the gas within 500 pc of the Sun 
are representative of the Galactic disk. Observations with the FUSE satellite should probe such 
more distant environments and hopefully answer some of these questions. 

The value of D /H for 7 2 Vel is larger than that usually considered typical for the Milky Way. 
This robust result establishes a new lower limit to the primordial D/H ratio. Within the framework 
of standard Big Bang Nucleosynthesis (Walker et al. 1991), the large value of D/H found toward 
7 2 Vel is equivalent to a cosmic baryon density of Hg/j 2 = 0.023 ± 0.002. This error simply reflects 
the uncertainty in the D/H ratio toward 7 2 Vel reported in this paper. We regard this value of 
f^/i 2 as an upper limit since no correction has been applied for the destruction of deuterium in 
stars. This upper limit on f^/i 2 is consistent with the preferred values of f^/i 2 derived from 
recent analyses of the BOOMERANG and MAXIMA cosmic microwave background measurements 
(e.g., Lange et al. 2000; Tegmark & Zaldarriaga 2000; Hu et al. 2000). However, any lowering of 
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this upper limit on f^/i 2 to correct for astration will lead to a marginal disagreement with simple 
inflation models (see Figure 4 in Tegmark & Zaldarriaga 2000) and requires adjustments of other 
cosmological parameters. Alternatively, the f^/i 2 upper limit from D/H may be taken as a prior 
assumption for the constraint of other cosmological parameters using the CMB data (e.g., see Hu 
et al. 2000). 
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